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PREFACE 


The  earth's  surface  may,  to  a  first  approximation, 
be  regarded  as  spherical  in  geometrical  terms  and,  in 
topological  terms,  it  may  be  considered  as  a  single 
continuous  closed  surface. 

In  physical  geography  the  sphericity  of  the  earth 
has  been  recognized  and  assigned  a  role.  Spatial  pre¬ 
diction  on  this  basis  has  been  so  very  successful  as  to 
have  been  commonplace  and  no  longer  regarded  to  be 
spatial  prediction  as  such.  Thus,  the  statements  about 
sunrise  and  sunset,  length  of  day,  seasonality,  celestial 
appearances  and  navigation,  and  the  like  which  depend 
in  part  upon  the  recognition  of  sphericity  of  the  earth's 
surface  are  regarded  as  certainties  and  the  predictive¬ 
ness  involved  is  lost  sight  of. 

The  topological  nature  of  the  surface  has  not  re¬ 
ceived  as  much  attention  and  the  requirements  that  this 
places  upon  the  theories  of  spatial  process  in  the  physical 
geography  of  global  patterns  of  circulation  have  not  been 
recognized,  to  the  disadvantage  of  those  theories. 

Imagine  now  two  spatially  continuous  distributions 
covering  the  surface  of  the  earth.  In  physical  geography 
wo  might  consider  temperature  and  pressure  distribution:*.. 

In  socio-economic  term:;  the  populat  ion  potential  field 
and  the  income  potential  field  serve  as  examples. 
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only  that  there  are  two  continuous  distributions,  A  and 
B,  over  the  earth's  surface,  there  will  always  be,  as 
Steinhaus  has  noted,  at  least  one  pair  of  antipodal  points 
having  both  the  samb  value  for  A  and  the  same  value  for 
B.  Through  time  the  spatial  distributions  involved  may 
change  and  the  given  antipodal  points  initially  so 
related  as  above  may  be  no  longer.  However,  another 
such  pair  will  then  exist  and  there  will  have  been  a 
shift  in  position.  This  is  the  consequence  of  the  theorem 
proposed  by  Ulam  and  proved  by  Borsuk  that  if  a  sphere 
is  folded  and  distorted  (but  not  torn)  so  as  to  be  made 
to  lie  flat,  there  is  necessarily  one  pair  of  antipodal 
points  from  the  sphere  which  come  to  lie  upon  each  other 
in  the  new  situation. 

The  proof  of  this  existence  theorem,  however,  pro¬ 
vides  no  means  for  finding  specific  locations  for  the 
equal-valued  antipodal  points  for  the  two  continuous 
distributions.  This  concerns  us. 

Theoretical  geography  is  a  science  of  earth  loca¬ 
tion  and  spatial  relations.  It  describes,  classifies, 
and  predicts  locations  in  the  spatial  sense.  Cartographies 
stands  to  geographical  science  as  graphics  docs  to  science 
generally.  Mapping  is  a  general  mathematical  concept 
in  the  theory  of  sets.  Cartography  is  the  geographical 
example  in  the  application  of  this  concept.  That  which 
is  ordinarily  called  a  map  by  geographers  and  laymen* is 
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technically  "a  graphical  image  of  a  mapping."  Whatever 
useful  roles  graphics  plays  in  science  generally  also 
can  be  claimed  for  cartography  with  relation  to  geograph¬ 
ical  science.  This  is  especially  true  now  that  geographers 
increasingly  employ  geometry  as  an-  appropriate  vehicle 
to  carry  their  discipline. 

There  is  already,  of  course,  an  accumulated  stock 
of  knowledge  and  experience  concerning  maps  as  stores  of 
spatially  ordered  information.  We  hope,  however,  to 
examine  the  expanded  roles  that  mapping  (especially 
computer-assisted  mapping)  seems  well-suited  to  play  in 
the  sciences  viewed  from  the  standpoint  of  theoretical 
geography  and  in  the  disciplines  employing  its  models 
for  decision-making  purposes. 

We  recognize  yet  another  role  for  maps.  In  the 
solution  of  certain  problems  for  which  the  mathematics , 
however  elegantly  stated,  is  intractable,  graphical 
solutions  (approximations)  are  possible.  This  is 
especially  true  with  regard  to  "existence"  theorems. 

There  arc  many  cases  in  which  the  graphical  solution  to 
a  spatial  problem  turns  out  to  be  a  map  in  the  full 
geographical  cense  of  the  term,  map.  Thus,  a  map  l_s 
the  solution  to  the  problem.  Computer-achieved  maps 
now  arc  possible  in  this  connection  whereas  previously 
computational  burdens  rendered  them  virtually  unattain¬ 
able. 
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GEOGRAPHY  AND  AN  EXISTENCE  THEOREM  --  A  Cartographic 
Computer  Solution  to  the  Localization  on  a  Sphere  of 
Sets  of  Equal-Valued  Antipodal  Points  for  Two  Contin¬ 
uous  Distributions  with  Practical  Applications  to  the 
Real  Earth 

One  can  discuss  the  existence  of  sets  of  antipodal 
points  on  spheres  without  concern  for  their  specific 
spatial  locations.  A  mathematician  working  with  a  certain 
existence  theorem,  for  instance,  may  be  primarily  inter¬ 
ested  in  whether  or  not  a  given  phenomenon  concerning 
the  pairing  of  antipodal  points  exists;  the  frequency 
of  occurence  and  location  may  be  beyond  the  scope  of 
the  proof.  Consider  the  following  example:  given  any 
two  continuous  functions,  f  and  g,  mapped  over  the  sur¬ 
face  of  a  sphere,  there  must  exist  at  least  one  set  of 
antipodal  points  at  which  function  f  is  equal  and  function 
g  is  equal.  The  proof  of  the  theorem  gives  ro  information 
about  how  many  other  such  pairs  of  points  exist  or  where, 
cn  the  surface  of  the  sphere,  any  of  such  sets  is  located. 
One  might  expect  that  geographers  and  others  would  be 
interested  in  questions  of  this  nature,  however.  This 
paper  is  an  attempt  to  answer  these  questions  by  devel¬ 
oping  a  procedure  to  find  spatial  solutions  to  the 
previously  stated  existence  theorem  with  the  aid  of 
and  existing  hYnAP  computer  mapping  program.  (»V.- 
paper  No.  3  of  the  Harvard  Paper.*»__li»  Theoretical  deog- 
ra£h y ,  "Implicit  map  Projections  in  Computer  Print-Pul;.," 
ami  paper  Ho.  li>,  "A  Two-pjr.cn::  jona.1  Interpol. at  if.-;-  I-Yr-cMe 
for  Computer  napping  of  Irregularly  .'“.paced  Pats,"  for  a 
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description  of  this  program.)  Stated  simply,  the  problem 
is  to  find  on  a  sphere  all  pairs  of  antipodal  points  6n 
the  sphere  at  which  the  values  for  some  continuous  function 
f  are  equal  and  the  values  for  some  continuous  function  g 
are  equal.  Of  course  the  values  for  each  of  these  functions 
may  differ  from  set  to  set  and  for  real  phenomena  the 
actual  values  depend  upon  the  units  chosen. 

The  first  step  in  rendering  the  problem  suitable  to 
a  computer  solution  was  to  change  the  problem  domain 
from  the  surface  of  a  sphere  to  a  plane.  There  are  any 
number  of  methods  of  projecting  a  sphere  onto  a  plane 
but  the  possible  choices  are  limited  depending  upon 
several  requirements  of  the  projection.  Most  importantly, 
the  projection  must  facilitate  the  comparison  of  function 
values  at  antipodal  points  on  the  sphere.  Also,  since 
the  continuity  of  the  functions  is  crucial,  the  projection 
should  make  a  minimum  of  cuts  in  the  surface  of  the 
sphere  and  those  should  be  of  a  regular  nature.  A  set  of 
two  hemispherical  projections  was  chosen  as  best  suited 
for  these  requirements.  If  the  sphere  is  split  into  two 
hemispheres  and  each  is  projected  onto  a  plane  lying 
parallel  to  the  cut,  a  rotation  of  one  of  the  hemispher¬ 
ical  projections  by  180“  will  bring  all  sets  of  antipodal 
points  into  alignment.  If,  however,  as  is  usually  the 
ca cc,  f  actions  are  made  from  antipodal  points  to 

tangent  pianc a,  on*  Of  these  be  inserted  before 

rotation  to  bring,  all  antipodal  points  into  alignment. 

Now  consider  what  the  locus  of  equal-valued  antipodal 
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points  night  look  like  on  the  hemispherical  projections. 

The  problem  is  further  simplified  by  assuming  the  globe 

♦ 

to  be  the  sum  of  many  great  circles  and  then  examining 
the  behavior  of  the  function  along  any  great  circle. 

The  conclusion  is  that  there  Is  at  least  one  set  of 
antipodal  points  where  the  function  has  equal  values  on 
every  great  circle.  The  proof  is  simple.  Since  the 
function  is  continuous  on  the  whole  sphere,  it  must  be 
continuous  on  the  great  circle.  Choose  a  diameter  of 
the  great  circle:  its  endpoints  are  antipodal  points. 
Rotate  the  diameter  along  itself  through  180°,  subtracting 
the  value  of  the  function  at  one  end  of  the  diameter  from 
the  value  at  the  oposite  end  (the  antipodal  point).  If 
the  initial  result  of  the  subtraction  be  some  number  A, 
when  the  diameter  has  completed  its  arc  the  result  of 
subtraction  would  be  -A,  Since  the  value  went  from  A 
to  -A  and,  because  the  result  of  subtraction  of  contin¬ 
uous  functions  is  continuous,  at-  some  set  of  points  the 
result  is  sero.  At  that  set  of  antipodal  points,  the 
function  has  the  same  value.  Thus  any  great  circle  has 
at  least  one  set  of  equal  antipodal  points. 

Superimpose  the  two  hemispherical  projections  of 
the  globe  so  that  each  point  on  one  hemisphere  is  aligned 
with  its  antipodal  point  on  the  other  hemisphere.  The 
result  is  a  circular  map,  each  of  whose  points  represents 
a  pair  of  antipodal  points,  Kvery  diameter  oi  this  circu¬ 
lar  map  is  o  greet  circle  of  tne  original  sphere.  (There 
are,  of  course,  other  great  circles  besides  these  diameters 
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The  characteristic  of  all  great  circles  on  these  circular 
maps  is  that  they  intersect  the  circumference  of  the  map 
at'  diametrically  opposite  points.)  It  was  previously 
shown  that  every  great  circle  (each  diameter  on  this  map) 
will  have  at  least  one  set  of  equal  antipodal  points.  If 
a  whole  series  of  diameters,  each  with  its  set  of  equal 
antipodal  points,  he  drawn  on  the  map,  a  reasonable  assump¬ 
tion,  by  virtue  of  the  continuity  of  the  function,  is 
that  the  sum  of  these  points  will  form  a  curve  of  some 
sort.  That  is,  in  the  limit  of  an  infinite  number  of 
diameter^  one  would  expect  the  locus  of  equal  sets  of 
antipodal  points  to  be  a  curve.  Note  that  where  the  curve 
intersects  the  circumference  (it  must,  since  the  circum¬ 
ference  is  a  great  circle)  it  roust  also  reappear  at  the 
diametrically  opposite  point.  This  is  because  -»ach  point 
on  the  circumference  net  only  represents  a  set  of  antipodal 
points,  but  is  the  antipodal  point  of  its  diametrically 
opposite  point  on  the  circumference  (i.e.,  the  circum¬ 
ference  duplicates  a  set  of  antipodal  points). 

Thus  the  locus  of  equal-valued  antipodal  point  sets 
should  be  a  closed  curve  (only  the  particular  projection 
used  makes  the  locus  appear  to  have  endpoints).  There 
can  also  be  a  single  isolated  set  e>f  equal-valued  antipod¬ 
al  points.  Those  events,  however,  represent  somewhat 
special  cases  and  will  be  discussed  later. 

One  might  also  ask  what  is  the  minimum  number  of 
such  closed  curves  and  what  is  tf.eir  local  Ion?  This 
minimum  locus  demands  that  every  conceive;  it  f.rc.i  circle 
intersect  it  at  lc:M  once.  Intuitively,  any  line* 


2.  The  locus  of  sets  of  antipodal 
points  with  equal  values  Is 
generated  point  t»y  point. 


connecting  two  opposite  points  on  the  circumference  of  . 
the  circular  map  will  Intersect,  at  least  once,  every 
great  circle  which  can  be  drawn.  This  follows  from  the 
previously  stated  observation  that  all  great  circles 
begin  and  end  at  pairs  of  diametrically  opposite  points 
on  the  circumference  of  a  circular  map.  Thus  the  sim¬ 
plest  locus  of  equal-valued  antipodal  points  might  be 
a  diameter  or  half  of  the  circumference  of  the  circular 
map,  or  more  generally,  the  projection  onto  the  circular 
map  of  any  great  circle  on  the  sphere.  This  is  a 
plausible  conclusion  since  every  great  circle  on  a  sphere 
intersects  every  other  great  circle  at  one  set  of  antip¬ 
odal  points. 

If  the  minimum  locus  of  equal-valued  antipodal 
points  for  one  function  is  any  line  connecting  two 
opposite  points  on  the  circumference  of  the  circular  map, 
the  same  must  be  true  for  any  other  function.  The  impos¬ 
sibility  of  connecting  each  of  two  sets  of  diametrically 
opposite  points  on  the  circumference  with  lines  which 
do  not  intersect  at  least  once  is  an  intuitive  affirma¬ 
tion  of  the  original  existence  theorem.  With  these 
expectations  as  to  the  form  of  the  solution,  the  solution 
procedure  can  now  be  developed. 

A  brute  force  solution  is  described  first.  For 
each  of  two  functions,  use  the  SYi’.AP  program  to  produce 
contour  maps  (on  hemispherical  projections  as  base  maps) 
of  both  hemispheres.  Taking  one  function  at  a  time, 
align  the  two  hemispherical  contour  maps  so  that  antipodal 
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points  are  aligned  and  mark  the  locations  where  the  same 
contour  levels  cross.  This  will  produce  not  the  line  as 
hypothesized  before,  but  a  string  of  polygons  whose  size 
will  depend  on  the  size  of  the  contour  levels.  (This 
will  be  discussed  in  more  detail  later.)  After  carrying 
out  the  same  procedure  for  the  second  function,  it  is  a 
simple  matter  to  compare  the  two  loci  of  equal-valued 
antipodal  points  to  see  where  they  intersect. 

The  first  step  is  to  produce  contour  maps  using  the 
SYMAP  program.  But  even  before  a  base  map  can  be  plotted 
a  particular  hemispherical  projection  must  be  chosen. 
Since  the  accuracy  of  the  SYMAP  interpolation  is  depen¬ 
dent  upon  the  area  over  which  the  interpolation  occurs, 
a  hemispherical  projection  which  preserves  the  area  of 
the  sphere  (although  this  will  mean  some  distortion  of 
linear  distances)  is  an  appropriate  choice. 

While  the  theorem  concerns  functions  mapped  on  any 
sphere,  the  examples  can  be  restricted  to  a  globe  of 
the  earth  without  loss  of  generality.  This  has  the 
advantage  of  providing  convenient  terminology,  such  as 
north  and  south,  and  a  grid  system  of  latitude  and 
longitude  for  locations  as  well.  Thus,  for  convon*  ice 
onl  •,  center  the  projections  on  the  North  and  South 
Poles;  the  circumference  of  the  previous  circular  maps 
is,  then,  the  equator. 

This  map  projection  is  called  a  polar  azimuthal 
equal  area  projection.  The  radius  of  a  latitude  line 
on  the  map  is  set  equal  to  the  length  of  the  chord  from 
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the  pole  to  the  same  latitude  line  on  the  globe.  Thus 
the  radius  of  the  projected  map  is  /2  times  the  radius 
of  the  globe  but  the:  area  of  the  map  is  equal  to  the 
area  of  one  hemisphere  of  the  globe.  The  greatest 
linear  distortion  will  occur  at  the  equator  where  the 
circumference  of  the  map  is  /?  times  the  circumference 
of  the  equator  on  the  globe. 

A  globe  of  4  inch  radius  was  chosen  from  which  to 
project  the  polar  equal  area  maps  (the  radius  of  the 
map  was  thus  5.66  inches).  These  numbers  were  chosen 
so  that  the  map  would  fit  on  one  panel  of  computer 
printout.  Choosing  the  number  of  data  point  locations  and 
spacing  was  a  less  rational  process.  The  data  points 
had  to  be  evenly  spaced  and  close  enough  together  to 
avoid  missing  phenomena  of  the  functions.  A  symmetrical 
arrangement  was  an  aid  to  plotting  data  point  locations 
and  facilitated  rotating  the  maps  so  that  the  antipodal 
points  were  aligned.  The  resultant  base  map  for  the  first 
series  of  maps  contained  6't  data  points  within  the  equator 
spaced  as  shown  on  Figure  4 . 

Functions  generated  with  random  numbers  were  used 
for  the  first  attempts  at  perfecting  a  method  to  solve 
the  problem  to  avoid  any  singularities  or  pecularities 
which  some  naturally  occur ing  function  might  have.  By 
assigning  random  numbers  from  1  to  99  to  every  data 
point  and  interpolating  with  SYY1A5*,  the  result  was  a 
contour  map  of  a  continuous  function,  with  one  Important 
exception,  i.e.,  at  the  equator.  If  random  numbers  were 
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Figure  3.  Projection  of  Polar  Azimuthal  Equal  Area 


Map  Prom  the  Globe 


assigned  to  data  points  in  vach  hemisphere  and  interpola¬ 
ted  out  to  the  equator,  the  function  values  at  the 
equator  of  both  hemispheres  were  generally  not  the  same. 
But,  the  functions  must  be  continuous  over  the  whole 
sphere  for  the  theorem  to  hold.  To  provide  continuity 
at  the  equator,  28  data  points  were  taken  from  one 
hemisphere  and  added  to  the  other  hemispherical  map  in 
the  region  beyond  the  equator.  For  example,  the  ring 
of  data  points  closest  to  the  equator  in  the  southern 
hemisphere  was  placed  beyond  the  equator  in  the 
appropriate  position  in  the  northern  hemispherical  map. 
The  computer  would,  then  interpolate  across  the  equator 
on  both  maps  and  provide*  the  same  values  along  the  equa¬ 
tor  in  both  hemispheres.  Subsequent  maps  confirmed 
that  values  along  the  equator  agreed  to  within  1  or  2 
character  locations  from  one  hemisphere  to  the  other. 

With  these  28  points  in  the  "overlap"  region,  there  was 
a  total  of  92  data  points  per  map. 

After  a  set  of  trial  maps  with  five  contour  levels 
confirmed  the  validity  of  the  solution  to  the  continuity 
problem  at  the  equator,  a  set  of  two  hemispherical  maps 
for  each  of  two  different  randomly  generated  functions 
was  produced  with  a  radius  of  11.32  inches  (two  panels 
of  computer  printout).  The  maps  were  done  in  10  levels 
with  white  contour  lines  between  the  levels.  These  had 
been  omitted  from  smaller  maps  because  they  constituted 
a  large  percentage  of  the  available  character  locations. 

The  procedure  for  finding  the  locus  of  equal-valued 
antipodal  points  is  as  follows.  For  ei.ch  function, 
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Generation  of  ’'Overlap"  Points 


Figure  6.  Random  Punctlon  Mapped  on  one  Hemisphere 
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place  a  sheet  of  clear  plastic  over  one  of  the  hemis¬ 
pherical  maps  and  trace  all  contour  levels  with  a  grease 
crayon,  labeling  the  contour  levels  with  appropriate 
values  (1-10).  Then  place  this  tracing  over  the  other 
hemispherical  map  so  that  antipodal  points  are  aligned 
and  mark  the  regions  where  contour  levels  of  the  same 
value  crdss.  The  result  for  one  of  the  randomly  genera¬ 
ted  functions  is  illustrated  below. 

As  noted  earlier,  the  locus  of  equal-vali  ?d  antip¬ 
odal  points  shows  not  as  a  series  of  closed  curves,  but 

% 

rather  as  a  string  of  polygons.  The  reason  for  this  is 
that  a  continuous  funclton  has  been  approximated  with  a 
step  function  and  all  points  in  the  same  contour  level 
are  taken  as  equal  in  value  when  in  fact  they  represent 
a  range  of  values.  As  the  number  of  contour  levels  in¬ 
creases,  the  range  of  each  contour  level  will  shrink. 
Because  the  functions  are  continuous,  in  the  limit  of 
an  infinite  number  of  step  functions,  each  of  vanishingly 
small  range,  one  would  expect  the  string  of  polygons  to 
shrink  to  a  series  of  closed  curves.  This  process  is 
illustrated  below  where  the  number  of  contour  levels  has 
been  doubled.  The  area  of  the  polygons  shrinks  approx¬ 
imately  by  a  factor  of  two.  Note  that  only  when  the  lines 
dividing  the  same  contour  levels  intersect  (point  A,  for 
example,  on  Figure  8a)  are  there  equal-valued  antipodal 
points.  Points  within  the  polygons  arc  sots  of  equal¬ 
valued  antipodal  points  only  to  within  the  value  range  of 
the  contour  level. 

The  strings  of  polygons  can  be  approximated  with  a 
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series  of  lines.  Although  the  precise  position  of  the, 
line  of  sets  of  equal-valued  antipodal  points  within  the 
polygons  is  often  uncertain,  there  is  a  series  of  points 
where  contour  lines  intersect  through  which  the  locus 
must  pass.  But,  there  is  still  considerable  uncertainty 
about  where  the  locus  of  equal -valued  points  is  located 
when  the  polygons  are  large  in  size  and  irregular  in 
shape.  (Figure  9b) 

Thus  it  is  with  uncertainty  that  the  strings  of 
polygons  are  reduced  to  a  series  of  lines  in  Figure  10. 

The  loci  of  sets  of  equal-valued  antipodal  points  for 
both  functions  are  drawn  on  the  same  map.  Within  the 
limits  of  the  aforementioned  uncertainty,  each  time  one 
of  the  loci  vanishes  at  the  equator  it  reappears  at  the 
opposite  point.  There  are  about  2^  points  of  intersection 
where  both  randomly  generated  functions  nave  equal  values 
at  sets  of  antipodal  points. 

This  method  gives  a  satisfactory,  though  not  ideal, 
approximation  to  the  problem.  By  using  the  computer 
maps  merely  as  stores  of  spatially  ordered  data,  problem 
solvirg  capabilities  have  not  been  utilized  to  the  fullest 
extent.  The  actual  manipulation  of  the  data  was  a  human 
effort  and  a  tedious  one  at  that.  As  a  refinement  of 
this  preliminary  technique,  the  computer  can  be  involved 
in  more  than  the  initial  display  of  the  data  —  it  can 
manipulate  the  data  in  some  way  relevant  to  the  solution 
of  the  problem  and  then  graphically  display  the  results  • 
of  the  manipulation. 
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Figure  9.  The  process  by  which  strings  of  oolyi'ons 
j  are  reduced  to  lines  Is  subject,  to  p;rcnt 
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Figure  10.  Loci  of  Sets  of  Rq uni  Valued  Antipodal  Points  for 
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The  first  step  is  to  eliminate  the  tedious  work 
involved  in  tracing  contour  levels  and  searching  for 
overlapping  contours,  a  process  that  accumulates  inac¬ 
curacies  with  multiple  tracings,  etc.,  and  is  susceptible 
to  errors  of  omission  on  large  complicated  maps.  The 
process  of  comparing  the  function  values  at  sets  of 
antipodal  points  is  analogous  to  subtracting  the  function 
value  at  each  point  from  the  value  at  its  antipodal  point. 
The  result  of  this  process  over  the  whole  hemispherical 
map  is  called  the  "difference  function",  a  continuous 
function.  Every  point  on  the  difference  function  map¬ 
ping  represents  two  points  on  a  globe,  a  point  on  one 
hemisphere  and  its  antipodal  point  on  the  other  hemis¬ 
phere.  When  the  difference  function  is  zero,  the  function, 
at  that  set  of  antipodal  points,  has  equal  values.  If, 
at  some  point,  the  difference  function  has  a  positive 
value  and,  at  a  nearby  point,  the  difference  is  negative, 
the  value  of  the  difference  furictio.  must  pass  through 
zero  on  any  path  Joining  the  two  points,  a  consequence 
of  the  continuity  of  the  function.  The  set  of. points 
where  the  difference  function  takes  on  the  value  of 
zero  is  the  locus  of  equal-valued  sets  of  antipodal  points. 

With  t!.;  use  of  the  computer,  the  thousands  cf 
subtractions  necessary  to  calculate  the  difference  function 
at  every  character  location  Is  not  a  difficult  feat.  A 
subroutine  can  bo  programmed  into  the  SYrlAP  program  which 
allows  the  computer  to  store  on  tape*  the  function  values  • 
at  every  character  location  for  two  different  functions 
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on  the  same  base  map  (actually  the  same  function  on  two 
different  hemispheres)  and  then  subtract  function  values 
at  each  character  location,  lie  computer  could  then  be 
instructed  to  print  symbolism  only  at  those  locations 
where  the  result  of  subtraction  was  zero,  or  since  a 
difference  of  exactly  zero  would  be  unlikely,  at  all 
locations  where  the  result  of  subtraction  was  between 
±6,  where  6  is  a  small  number  which  would  vary  with 
each  function. 

As  first  approximation  to  subtraction  at  every  char¬ 
acter  location,  several  difference  functions  were  generated 
by  computing  the  value  of  the  difference  function  only 
at  data  points  (by  subtracting  the  value  at  a  data 
point  in  the  northern  hemisphere  from  the  value  at  the 
antipodal  point  in  the  southern  hemisphere)  and  then 
interpolating  over  the  whole  map  to  produce  values  for 
the  difference  function  at  every  character  location. 

This  was  considerably  simpler  and  less  costly  than  writing 
a  subroutine,  debugging  it  and  then  subtracting  at  about 
8000  character  locations. 

Difference  functions  were  mapped  in  10  levels  with 
no  white  spacing  for  contour  lines  for  each  of  the  two 
randomly  generated  functions,  had  there  been  many  data 
points  or  large  numbers  to  subtract,  the  values  of  the 
function  at  corresponding  antipodal  points  in  each 
hemisphere  could  have  bee.  entered  in  a  data  bank  and 
the  subtraction  performed  by  the  computer.  Since  there 
were  only  9?  subtractions  of  two  digit  numbers,  it  was 


Ftp.urc  11.  The  aero  contour  for  the  difference  function  i 
marked  between  levels  l>  and  6.  The  equate  r  lr» 
also  indicated. 


just  as  simple  to  do  the  subtractions  by  hand  and  use 
the  results  as  data  point  values  of  the  difference 
function  mapping. 

The  printout  of  the  0  level  for  one  of  the  difference 
functions  is  shown  in  Figure  11.  For  this  function,  the 
0  level  corresponded  to  a  contour  line  between  levels  5 
and  6  and  thus  could  be  easily  plotted.  For  the  second 
function,  the  0  level  contour  fell  in  the  middle  of  a 
contour  level  so  there  was  a  certain  amount  of  guesswork 
as  to  where  it  actually  was.  The  resultant  locus  of 
equal-valued  sets  of  antipodal  points  for  each  of  the 
two  functions  was  in  general  agreement  with  the  previously 
plotted  results. 

Having  now  established  two  preliminary  methods  for 
finding  the  location  of  the  points  that  the  existence 
theorem  says  must  exist,  I  decided  to  do  further  work 
with  real  functions  as  opposed  to  the  functions  generated 
with  random  numbers.  Two  completely  unrelated  functions 
could  have  been  used  and  the  theorem  would,  of  course  have 
held,  but  if  for  real  phenomena  related  functions  were  used 
there  could  have  been  some  significance  to  the  spatial  loca¬ 
tions  of  the  sets  of  antipodal  points  at  which  each  function 
had  equal  values.  The  choice  of  functions  was  governed 
mainly  by  what  was  available.  1  desired  two  related 
functions  mapped  over  the  whole  earth  with  a  minimum  of 
10  contour  levels  and  preferably  displayed  on  a  polar 
equal  area  projection.  Two  maps  in  the  series  "World 
Maps  of  Climatology"  by  H.K.  Lands'oerg,  H.  Lippman,  K.H. 
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Paffem,  and  C.  Troll,  1963,  met  the  first  two  requirements 
but  not  the  last.  The  two  functions  were  total  sunshine 
received  by  the  earth  in  hours  per  year  (hereafter  called 
the  Sunshine  function)  and  average  annual  energy  from 
sky  radiation  in  kcal/sq.  cm/year  (the  Energy  function). 

Were  it  not  for  atmospheric  factors,  etc.,  both  functions 
would  have  depended  on  latitude  only  and  the  functions  would 
have  had  equal  values  at  every  set  of  antipodal  points. 
Because  of  atmospheric  and  other  disturbances,  the  maps 
showed  great  deviations  from  a  latitude  dependent  phenomena. 
It  might  be  revealing  to  note  the  sets  of  antipodal 
points  at  which  the  factors  which  made  one  function  have 
equal  values  also  made  the  second  equal-valued.  This  was 
just  a  restatement  of  the  original  theorem. 

Since  the  Energy  and  Sunshine  function  were  not 
available  on  hemispherical  maps,  the  contours  were  redrawn 
on  a  polar  projection  to  facilitate  finding  the  function 
values  at  data  points.  This  was  not  really  necessary  — 
the  location  of  each  data  point  on  the  polar  projection 
could  have  been  determined  in  terms  of  latitude  and 
longitude  and  then  the  function  value  at  that  point  on 
the  original  map  could  have  been  found.  One  could  also  have 
derived  equations  to  convert  from  latitude  and  longitude 
on  a  sphere  to  SYMAP  coordinated  on  the  base  map.  The 
equations  arc: 
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for  a  polar  azimuthal  equal  area  projection  from  a  globe 
of  radius  R  where  4>  =  latitude,  0  «*  longitude  and  (x0  ,y0 ) 
are  coordinates  of  the  center.  (0  is  measured  counter¬ 
clockwise  from  the  top  of  the  map.)  But  this  would  not  be 
really  useful  unless  one  were  willing  continually  to  change 
data  point  locations. 

The  loci  of  equal-valued  sets  of  antipodal  points 
for  the  Energy  and  Sunshine  functions  were  determined 
first  by  the  laborious  hand  tracing  technique  and  then  by 
using  the  difference  functions.  As  before,  while  there 
was  some  uncertainty  in  the  exact  position  of  the  equal- 
value  loci  with  the  hand  tracing  method,  every  time  the 
locus  intersected  the  equator  it  reappeared  on  the  opposite 
side  as  it  should  have.  The  difference  function  approach 
lacked  this  virtue  although  when  the  zero  contour  fell 
between  two  contour  levels  there  was  no  guesswork  as  to 
the  location  of  equal-valued  sets  of  antipodal  points. 

It  appeared  that  the  difference  function  missed  some  of 
the  0  contours  because  the  data  point  spacing  was  too  wide. 

As  a  refinement  of  this  technique,  reconsideration  of 
the  possibility  of  subtracting  at  every  character  location 
was  in  order.  An  examination  of  the  interpolation  process  elim¬ 
inated  this  process  for  the  following  reasons. 

Subtraction  at  every  character  location  will  produce 
a  significant  improvement  in  the  accuracy  only  if  the 
function  values  to  be  subtracted  have  been  computed  to 
great  accuracy  at  every  character  location.  The  SYilAF  program 


The  shaded  area  is  the  locus  of  sets  of  antipodal 
points  with  equal  function  values.  The  function  i 
the  total  annual  cnerrv  Incident,  upom  the  earth's 
surface . 


The  shaded  area  Is  locus  of  sets  of  antioodal 
points  with  equal  function  values.  The  function 
is  the  annual  sunshine  duration  at  the  earth’s 
surface. 
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Fir.ure  12d. 


The  zero  contour  for  the  sunshine  difference 
function  is  drawn. 
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instructs  the  computer  to  do  sophisticated  Interpolation  at 
every  third  character  location  in  the  horizontal  direction 
and  every  other  character  location  in  the  vertical  direction 
and  then  uses  simple  linear  interpolation  to  fill  in  the 
gaps.  Thus  subtraction  at  every  character  location  is 
not  necessary  —  the  result  would  be  virtually  indistin¬ 
guishable  from  placing  data  points  at  every  third  location 
and  then  interpolating  between.  In  fact,  the  resulting 
map  would  not  differ  significantly  if  even  fewer  data 
points  were  used  and  this  would  save  the  trouble  and 
expense  of  writing  an'*  debugging  a  subroutine. 

Thus,  to  refine  the  technique  further,  subtraction 
at  every  character  location  was  bypassed,  although  that 
would  certainly  have  worked  and  have  been  an  improvement. 
Instead,  a  new  base  map  was  layed  out  with  almost  three 
times  as  many  data  points  (a  total  of  201  including  overlap 
points)  and  an  average  spacing  between  adjacent  data 
points  of  about  nine  character  locations.  The  equator 
was  the  same  size  as  before  and  was  approximated  by  a 
32  sided  polygon.  The  series  of  overlap  points  beyond 
the  equator  remained  to  insure  continuity  at  the  equator 
and  a  series  of  data  points  lay  on  the  equator  itself 
(on  the  previous  base  map  there  were  no  data  points  on 
the  equator).  These  improvements  were  designed  to  overcome 
the  omissions  of  previous  difference  mappings.  Finally, 
on  this  series  of  maps  the  0°  longitude  line  was  rotated 
by  l0O°  to  the  top  of  the  map  and  the  symbolism  ended 
at  the  equator,  not  at  the  overlap  points  as  on  the 
previous  maps. 
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Figure  13.  Data  Point  Location*  on  one  Quadrant  of  Revised 
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Figure  1<I .  Revised  Rase  Map  for  Mapping  Difference 
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To  eliminate  the  guesswork  involved  in  plotting  the 
0  contour  when  it  falls  in  the  middle  of  a  contour  level, 
several  of  the  electives  available  with  SYMAP  are  utilized. 
These  electives  are  used  to  make  the  computer  actually 
print  out  the  zero  contour  and  thus  eliminate  another 
human  task  in  the  solution  procedure.  This  is  done  by 
specifying  a  map  of  the  difference  function  in  two  levels 
only,  with  the  zero  contour  line,  or  locus  of  equal¬ 
valued  sets  of  antipodal  points,  between  them.  Appro¬ 
priate  symbolism  might  be  +  signs  for  all  regions  where 
the  difference  functions  is  greater  than  zero  and  -  signs 
for  all  regions  less  than  zero.  The  darkest  symbol  is 
specified  for  the  contour  line  to  make  it  stand  out. 

These  changes  would  produce  a  highly  readable  difference 
map  and  locate  the  0  contour  to  within  the  accuracy  of 
the  original  data.  But,  even  more  information  can  be 
obtained  from  the  same  map.  Since  there  will  be  some 
error  involved  in  estimating  the  function  values  on  the 
original  contour  map,  the  spatial  location  of  the  0 
contour  may  not  be  located  precisely.  An  error  in 
estimating  the  original  function  value  would  show  up  as 
a  small  displacement  of  the  0  contour  on  the  difference 
map.  An  estimate  may  be  made  of  just  how  far  the  0 
contour  might  stray  from  its  calculated  position  by 
putting  limits  on  the  values  permitted  in  each  contour 
level.  If  an  a  ’age  error  In  reading  an  Energy  value 
Is  taken  as  +  kcal/cm  /year  then  the  propagated  error 
in  the  difference  function  is  +  7  kcai/cm  /year.  Thus, 
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lf  the  range  of  the  upper  contour  level  is  limited  to  +7 
and  the  lower  contour  level  to  -7,  the  appropriate  symbolism 
will  be  printed  only  where  the  value  of  the  difference 
function  is  within  those  ranges.  All  other  areas  will  be 
blank.  The  computer  output  would  consist  of  black  Q 
contour  lines  surrounded  by  bands  of  symbolism  of  various 
widths  within  which  the  0  contour  might  actually  be  located, 
if  there  were  some  error  in  estimating  function  values. 

Using  certain  mapping  electives  to  produce  the 
desired  results,  refined  difference  mappings  were  made  of 
the  Energy  and  Sunshine  functions  (average  propagated 
error  for  the  Sunshine  difference  mapping  was  estimated 
at  +  140  hrs/yr).  The  results  were  excellent.  In  most 
regions  of  the  map  the  contour  levels  were  wide  enough  to 
permit  the  black  zero  contour  to  be  printed.  In  other 
regions  where  the  contour  levels  were  only  one  character 
space  wide,  (A,  Figure  ibb),  the  ’else  location  of  the 
zero  level  was  known  (between  adjacent  +  and  -  signs) 
even  though  the  black  contour  was  omitted.  Only  where 
the  difference  function  was  changing  so  quickly  that 
scattered  symbolism  appeared  (13,  Figure  15b)  was  there 
some  uncertainty  about  the  exact  position  of  the  zero 
contour,  but  even  this  uncertainty  was  small.  Each  time 
the  zero  contour  intersected  the  equator  it  reappeared  at 
the  diametrically  opposite  point.  This  former  problem 
area  was  finally  eliminated.  Where  previously  the  exact 
behavior  of  the  zero  contour  had  been  in  doubt  It  was 
precisely  defined. 


F.nerpty  Function 


The  zero  contour  displays  many  of  the  properties  of 
the  intersection  of  two  surfaces.  If  the  difference  . 
function  is  envisioned  as  a  three-dimensional  surface, 
the  zero  contour  may  be  obtained  by  intersecting  the 
surface  with  a  level  surface  of  zero  value.  Although 
all  the  zero  contours  of  both  functions  form  a  series 
of  closed  curves  (ignoring  discontinuities  at  the  equator) 
the  possibility  of  Isolated  points  or  line  segments  at 
the  zero  level  must  not  be  overlooked.  Under  what 
conditions,  then,  do  two  continuous  surfaces  intersect 
in  a  point  and  in  a  line' segment? 

Point  intersection  may  be  interpreted  as  the  limiting 
stage  of  intersection  between  a  level  surface  and  a  peak. 
When  the  surface  intersects  the  base  of  the  peak,  the 
intersection  is  a  closed  curve.  As  the  level  surface  is 
moved  upward,  however,  the  closed  curve  shrinks  in  size 
to  a  point  and  then  vanishes  (Figure  l6a).  It  would  appear 
from  this  example  that  the  condition  for  a  single  point 
of  intersection  of  the  two  function  surfaces  (i.e.,  the 
locus  of  equal-valued  sets  of  antipodal  points  is  just 
one  set  of  points)  is  that  the  tangent  level  surfaces 
to  each  function  at  that  particular  point  be  Identical. 

The  case  of  the  line  segment  can  be  evaluated  in  a  similar 
manner.  In  terms  of  intersecting  surfaces,  consider,  for 
example,  a  horizontal  surface  intersecting  a  ridge  with 
a  horizontal  ridgeline.  As  the  level  surface  moves  up 
the  ridge,  the  intersection  shrinks  from  an  Initially 
elongated  clo-eu  curve  to  n  line  segment  when  the  surface 

is  Just  tangent  at  the  ridgeline.  As  previously  mentioned, 
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neither  the  single  point  nor  line  segment  loci  of  equal- 
valued  sets  of  antipodal  points  were  observed  on  these 
two  difference  function  mappings. 

The  particular  printout  format  of  the  difference 
functions  has  yet  another  potential  use.  It  provides 
information  which  aids  in  answering  the  following  question: 
What  happens  to  the  spatial  solutions  of  the  existence 
theorem  as  one  or  both  of  the  functions  is  varied  by 
small  amounts?  As  the  first  step  in  answering  this 
question,  the  spatial  variation  of  the  zero  contour  of 
one  of  the  difference  functions  under  small  perturbations 
of  the  function  is  examined. 

Only  the  case  of  small  perturbations  of  the  function 
in  a  localized  region  is  reviewed  because  larger  variations 
would  result  in  changes  toe  widespread  and  complex  to 
analyse.  They  would  present,  essentially,  a  completely 
new  function  and  could  be  treated  as  such,  by  mapping  the 
difference  function,  etc.  The  distinction  between  large 
and  small  variations  of  the  function  is  vague.  It  will 
be  convenient  to  consider  small  perturbations  as  those  on 
the  order  of  the  range  of  tne  contour  levels  of  the 
difference  function.  This  implies,  of  course,  that  what  • 
is  considered  small  for  one  function  may  be  large  for 
another,  s’nco  the  contour  level  ranges  will  vary  consid¬ 
erably  among  difference  functions. 

The  computer  printout  of  the  zero  contours  for  a 
difference  fuuc t lor.  yields  much  information  on  this  prob¬ 
lem.  ihe  total  width  (measured  perpendicular  to  ft 
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tangent  to  the  contour  at  that  point)  of  the  two  contour 

levels  at  any  particular  point  on  the  zero  contour  gives 

a  rough  indication  of  the  slope  of  the  difference  function 

at  that  set  of  antipodal  points.  First,  observe  that  the 

change  in  the  difference  function  (Af)  over  the  two  contour 

levels  is  just  twice  the  propagated  average  error  of  the 

original  function  (call  this  2c).  An  approximate 

corresponding  change  in  the  domain  of  the  function  (Ax) 

is  determined  by  actual  measurement  of  the  width  of  the 

two  contour  levels.  (Note:  the  scale  of  the  map  is 

% 

constant  in  a  circumferential  direction  but  varies  in 
the  radial  direction.)  Thus  the  gradient  of  the  difference 
function  can  be  approximated  at  any  point  on  the  zero 
contour  by 

To  avoid  complications  which  arise  from  the  varying 
scale  of  the  difference  map,  the  gradient  of  the  difference 
function  will  be  discussed  in  qualitative  terms  only. 

Where  the  gradient  is  small,  (that  is,  where  the  width 
of  the  two  contour  levels  is  large  (A,  Figure  17)  a 
small  change  in  the  value  of  the  difference  function  at 
some  point  on  the  zero  contour  would  result  <n  a  relatively 
large  change  in  the  spatial  location  of  the  zero  contour. 
For  instance,  if  the  function  is  varied  so  that  the 
difference  function  changes  by  "t"  at  some  point  on  the 
zero  contour,  one  would  expect  the  zero  contour  to  shift 
to  one  side  by  about  the  width  of  the  contour  level. 

The  direction  of  the  shift  would  clipen**  upon  wh».  th*  r  the 
difference  function  increased  by  c  or  d-»  c reaped  by  r. 
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Kipure  17.  hispl aceinont  of  the  zero  contour  of  the?  difference 
function  for  a  srrail  function  chan  re  at  so  no  set.  of 
antipodal  points  varies  with  the  inverse  o°  the 
gradient  at  the  points  In  question. 
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Similarly,  in  a  region  where  the  gradient  is  steep 
(B,  Figure  17)  the  difference  function  is  changing 
rapidly.  If  the  difference  function  is  changed  by  e  at 
some  point  on  the  zei  o  contour  in  this  region,  the  zero 
contour  will  experience  a  much  smaller  displacement  (one 
character  location  or  less). 

Thus  the  displacement  of  the  zero  contour  of  the 
difference  function  as  a  result  of  small,  localized 
function  variations  can  be  estimated  by  examining  the 
printout  as  described  above.  These  conclusions  hold  for 
changes. in  the  difference  function  resulting  from  small 
spatial  shifts  of  the  original  function  as  well  as  small 
variations  in  the  original  function  values.  The  continuity 
of  the  function  assures  the  equivalence  of  these  two 
cases. 

Only  small  changes  in  the  difference  function 
have  been  discussed  because  the  primary  concern  here  is 
with  what  happens  to  the  zero  contour  In  a  small  region 
about  the  point  in  question.  A  large  change  in  the  dif¬ 
ference  i unction  at  a  given  point  would  produce  significant 
changes  in  the  locality  under  consideration  and  the 
approximate  techniques  discussed  previously  would  not 
give  meaningful  results.  But,  because  the  distinction 
between  large  and  small  changes  is  necessarily  vague 
and  arbitrary,  it  may  at  some  time  be  necessary  to 
approximate  the  displacement  of  the  zero  contour  for 
difference  function  changes  as  large  as  2c  or  3c.  This 
can  be  accomplished  with  only  slight  revisions  in  the 
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previous  techniques.  By  changing  the  electives  which  control 
the  number  of  contour  levels  and  their  ranges,  a  difference 
map  could  easily  be  produced  with  four  levels,  for  instance, 
which  would  cover  twice  the  range  of  previous  maps.  To  use 
the  Energy  function  as  an  example,  the  levels  would  range 
from  14-7,  7-0,  0-(-7),  and  (-7M-14).  The  displacement 
of  the  zero  contour  could  then  be  estimated  in  the  locality 
of  a  point  where  the  Energy  difference  function  changed  by 
as  much  as  14  keal/cm  /year. 

If  the  approximate  displacements  of  the  zero  contours 
for  small  changes  in  both  functions  are  known  in  some  region, 
the  behavior  of  an  intersection  of  the  zero  contours,  i.e., 
the  pair  of  antipodal  points  at  which  each  function  has  equal 
values,  can  then  be  inferred.  The  gradient  of  the  difference 
function  at  any  point  on  the  zero  contour  line  is  perpendicular 
to  the  tangent  to  the  contour  at  that  point .  Small  changes 
in  the  difference  function  will  displace  the  zero  contour 
roughly  parallel  to  its  former  position  in  the  direction 
of  either  plus  or  minus  the  gradient  at  any  particular  point 
on  the  contour.  Thus,  to  approximate  the  motion  of  the 
point  of  intersection  of  two  zero  contours,  take  the  vector 
sum  of  the  Inverse  of  the  gradient  of  each  contour  (^~)  at 
that  point.  (The  direction  is  determined  by  the  sign  of  the 
change.)  The  resultant  vector  indicates  the  instantaneous 
direction  and  magnitude  of  motion  of  the  point  of  inter¬ 
section,  An  estimation  of  its  displacement  may  be  made  by 
plotting  the  new  locations  oi  both  zev,o  contours  and  then 
noting  the  new  position  of  their  point  of  Intersection, 
assuming  they  still  intersect. 


Points  of  intersection  of  zero  contours  may  vanish 
when  the  difference  functions  are  varied  slightly  but 
the  theorem  requires  the  existence  of  at  least  one  point 
of  intersection  regardless  of  how  greatly  the  functions 
change,  as  long  as  the  changes  are  continuous.  The 
significance  of  the  particular  spatial  locations  of  the 
sets  of  antipodal  points  at  which  both  of  the  functions 
have  equal  Values  is  a  complex  problem.  There  is  an 
unclear  dependence  on  the  degree  to  which  the  two  functions 
are  related.  The  Energy  and  Sunshine  functions,  for 
instance,  are  in  some  ways  related  (one  would  expect  the 
Energy  function  to  he  dependent  on  the  Sunshine  function) 
and  thi ^  is  borne  out  by  a  superposition  of  the  zero 
contours  of  both  difference  functions.  The  zero  contours 
of  these  two  difference  functions  intersect  over  twenty 
times  but  this  in  itself  is  not  significant  because  the 
randomly  generated  functions  previously  developed  had 
over  twenty  Intersections  as  well.  What  is  more  signif¬ 
icant  is  the  degree  to  which  the  zero  contours  of  one 
function  seem  to  correspond  to  those  of  the  other.  In 
many  regions,  the  contours  are  nearly  parallel  but 
slightly  displaced  —  in  others  they  practically  coincide. 
These  results  were  not  totally  unexpected. 

As  previously  mentioned.  Sunshine  and  Energy  are 
basically  latitude  dependent  phenomena.  Wore  i*  not  for 
the  presence  of  atmospheric  phenomena,  which  in  turn  are 
in  many  ways  dependent  on  other  physical  factors,  the  zero 
contours  of  both  difference  functions  would  have  wavered 


Fir.urc  lOa.  Loci  of  Sets  of  Kauai  Valued  Antipodal  Points  for 

Kncrr.y  Function  (Thick  Lines)  and  Sunshine  Function 


(Thin  Liner- ) 

note.- :  The  1  ntorr'’Ct  5 on:-  arc  circled. 


I'M r upo  lHh.  Kt.ro  Contours  of  Knorr.v  PI  f fVrenoe  Pune 


Line;! )  and  PunsMne  plfTei 


ICO  VUMct.  U'M 


-48- 

the  whole  sphere.  Whers  the  zero  contours  do  in  fact  exist 
may  indicate  that  the  same  atmospheric  conditions  and  their 
preconditions  exist  at  those  sets  of  antipodal  points. 

To  follow  up  this  hunch,  a  difference  function  was 
mapped  for  a  continentality  potential,  computed  by  Donald 
Shepard.  He  divided  the  earth  into  trapezia  of  5°  x  5° 
and  assigned  a  value  to  each:  100  if  the  trapezoid  con¬ 
tained  mainly  land,  0  if  it  was  mostly  water.  A  potential 
value  was  then  calculated  for  each  trapezoid  by  weighting 
the  value  of  every  other  trapezoid  by  the  inverse  of  the 
distance-squared.  A  contour  line  of  value  about  70  traces 
out  recognizable  continents.  Thus,  the  continentality 
potential  closely  resembles  a  map  of  the  earth,  accen¬ 
tuating  lines  of  demarcation  between  water  and  land 
masses. 

Since  the  continentality  potential  is  a  continuous 
function  over  the  surface  of  a  sphere  the  zero  contours 
of  its  difference  function  must  meet  the  criteria  pre¬ 
viously  established.  Examination  of  the  difference 
function  shows  that  it  does.  This  difference  function  is 
Interesting  because  it  has  very  little  more  than  what  was 
determined  as  the  minimum  locus  of  sets  of  antipodal 
points  with  equal  value;  that  is,  a  curve  connecting 
two  opposite  points  on  the  equator  o.  the  hemispherical 
map.  There  are  only  two  small  additional  closed  curves. 

Figures  20a  and  20b  show  the  zero  contours  of  the 
continentality  potential  difference  function  plotted  (on 
the  same  hemispherical  map  a?  the  Energy  and  Sunshine 


Fir.ure  ?0b.  Zero  Contours  of  Sunshine  n.l f  femnee  Function  (Thin 
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dlf.ference  functions.  Although  there  are  few  points 
of  intersection,  large  segments  of  the  continentality 
zero  contours  closely  approximate  the  zero  contours 
of  both  the  Energy  and  Sunshine  difference  functions. 

This  would  seem  to  imply  that  there  is  some  correlation 
between  the  factors  that  make  the  Energy  and  Sunshine 
functions  have  equal  values  at  antipodal  points  and 
the  degree  of  continentality  and  marine  influence  at 
these  sets  of  antipodal  points. 

Furthermore,  if  the  geographical  locations  of 
the  zero  contours  of  these  three  difference  functions 
are  examined,  there  appears  to  be  more  than  coinciden¬ 
tal  agreement  with  the  land-water  boundaries  on  the 
globe.  This  suggests  that  those  atmospheric  phenomena 
which  make  the  Energy  and  Sunshine  functions  equal  at 
sets  of  antipodal  points  are  the  result  of  physical 
conditions  related  to  the  influence  of  land  and  water 
masses  at  the  boundaries  between  them.  Couched  in  such 
general  terms,  this  proposition  may  be  a  statement  of 
the  obvious,  but  the  insights  developed  by  this  par¬ 
ticular  method  of  analysis  might  produce  valuable 
results  with  more  intensive  study. 

It  is  difficult,  if  not  Impossible,  to  assign  great 
significance-  to  the  geographical  location  of  the  inter¬ 
section  of  zero  contours  of  two  functions  if  the  functions 
are  totally  unrelated.  The  loeui*  of  intersections  would 
be  little  more  than  a  mathematical  curiosity.  An  attempt 
to  fabricate  causal  relations  among  physical  phenomena 
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to  explain  the  intersections  of  totally  unrelated  functions 

might  te  an  exercise  in  futility  although,  of  course,' 

the  results  might  be  suggestive  of  patterns  of  relationships, 

hitherto  ignored.  Caution  should  be  exercised  when  trying 

to  explain  the  intersection  of  zero  contours  of  functions 

which  are  related  in  an  unclear  way.  The  mathematical 

necessity  of  at  least  one  point  of  intersection  limits 

the  significance  one  should  attach  to  any  explanation  which 

ascribes  to  other  factors,  such  as  geological,  climatic, 

etc.,  the  responsibility  for  the  points  of  intersection. 

% 

These  other  factors  do  play  a  role,  perhaps,  in  determin¬ 
ing  the  number  of  intersections  and  .aeir  locations,  but 
this  depends  to  a  large  extent  upon  the  nature  of  the 
particular  functions  under  consideration.  On  the  other 
hand  theories  ignoring  the  necessity  of  these  intersections 
are  similarly  suspect. 

The  uncertainties  and  inaccuracies  resulting  from 
repeated  manipulation  of  the  original  data  have  been 
continually  stressed.  In  Judging  the  reliability  of 
the  final  results,  one  must  also  consider  the  Inherent 
limitations  in  the  use  of  the  mapping  program.  It  is 
critically  important  that  the  spacing  of  the  data  points 
be  finer  than  the  scale  of  the  phenomena  which  is  being 
mapped.  This  problem  plagued  the  first  attempt  to  con¬ 
struct  a  difference  function.  The  data  point  spacing  was 
so  wide  that  two  adjacent  data  points  might  be  in  positive 
regions,  for  instance,  while  a  negative  region  between 
them  was  missed,  resulting  in  the  omission  of  a  zero 


contour  line.  This  situation  was  rectified  by  the  addition 
of.  more  d&tapoints.  There  could  conceivably  be  a  function 
for  which  the  addition  of  more  data  points  would  still 
not  solve  this  problem.  A  detailed  function  showing 
population  potential  of  the  world  would  have  "spikes"  at 
most  city  locations.  But  even  a  reasonably  sized  SYMAP 
contour  map  with  data  points  at  every  character  location 
would  probably  miss  a  great  many  of  the  city  spikes.  If 
the  intention  were  to  show  the  general  character  of  the 
potential  then  the  map  would  have  sufficed,  but  if  it 
were  important  to  show  the  precise  function  value  at 
city  locations,  a  considerably  larger  map  would  have  had 
to  be  used. 

The  solution  procedure  advocated  in  the  preceding 
pages  for  finding  spatial  solutions  for  the  original 
existence  theorem  is  not  necessarily  the  best  solution  or 
a  final  solution  in  any  sense.  It  is,  however,  a  relatively 
quick,  accurate  and  inexpensive  computer  solution  to  the 
problem. 

To  review  briefly,  a  base  map  was  constructed  from 
polar  azimuthal  equal  area  projections  of  hemispheres  of 
a  globe  such  that  each  point  on  the  map  represented  a 
pair  of  points  on  a  globe  —  a  given  point  in  one 
hemisphere  and  its  antipodal  point  in  the  other  hemis¬ 
phere.  (Antipodal  points  were  aligned  by  a  simple  rota¬ 
tion  of  one  of  the  hemispheres.)  The  base  map  contained 
a  large  number  of  closely  spaced  data  points  and  a  set 
of  points  beyond  the  equator  to  insure  the  continuity  of 


the  function  across  the  equator.  The  "difference  function 
of  the  function  under  consideration  is  plotted  on  this 
base  map.  The  value  of  the  difference  function  at  a 
given  set  of  antipodal  points  is  computed  by  subtracting 
the  value  of  the  function  at  the  appropriate  point  in 
the  northern  hemisphere  from  the  value  at  its  antipodal 
point  in  ♦he  southern  hemisphere.  The  SYMAP  program  is 
then  used  to  produce  a  contour  map  of  the  difference 
function.  The  zero  values  of  the  difference  function 
are  of  prime  importance  since  at  these  sets  of  antipodal 
points,  the  original  function  will  have  equal  values. 

Using  the  appropriate  electives  in  the  SYMAP  program,  the 
computer  is  instructed  to  print  out  a  black  contour  line 
at  the  value  zero  and  one  contour  level  on  each  side  of 
the  zero  contour.  The  range  of  the  contour  levels  is 
determined  by  estimating  the  average  error  in  assigning 
function  values  at  data  points.  Where  the  difference 
function  changes  so  rapidly  that  there  are  insufficient 
character  locations  to  print  the  black  zero  contour, 
it  must  lie  between  adjacent  positive  and  negative  valued 
points.  The  two  contour  levels  not  only  provide  an 
estlmace  of  the  possible  inaccuracies  in  the  location  of  * 
the  zero  contour  but  also  give  an  indication  of  the 
gradient  of  the  difference  function.  This  allows  evalua¬ 
tion  of  the  displacement  of  the  zero  contour  for  varia¬ 
tions  in  the  original  function  whi  esult  in  small 
changes  in  the  difference  function. 

To  find  where  both  functions  each  have  equal  values 
at  sets  of  antipodal  points,  the  zero  contours  of  the 
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dlfference  functions  of  the  twc  original  functions  were 
drawn  on  the  same  map  and  the  points  of  intersection  were 
marked.  These  are  the  desired  points.  The  displacements 
of  the  points  of  intersection  resulting  from  variations 
in  both  original  functions  can  be  estimated  by  noting  the 
displacement  in  both  zero  contours  and  finding  the  new 
point  of  intersection  (assuming  they  still  intersect 
after  having  been  displaced). 

Comparison  of  results  obtained  by  this  method  and 
by  the  tedious  process  of  hand-tracing  and  overlapping 
the  original  contours  showed  that  none  of  the  zero  contours 
had  been  omitted  but  that  many  ambiguities  and  uncertain¬ 
ties  inherent  in  the  tracing  and  overlapping  method  had 
been  eliminated.  Improvements  in  this  procedure  can 
still  be  made.  If  the  accuracy  of  the  data  be  sufficiently 
good,  function  values  could  be  subtracted  at  every  char¬ 
acter  location  rather  than  only  at  data  points.  As 
stated  previously,  however,  because  of  the  nature  of 
the  interpolation  process,  this  method  would  not  produce 
significant  improvements.  The  use  of  more  sophisticated 
computing  equipment  which  produces  line  drawings  of  three- 
dimensional  surfaces  on  paper  or  cathode  ray  tubes  might  • 
vastly  improve  comprehension.  Superposition  of  function 
surfaces  or  intersection  of  the  difference  function  surface 
with  a  zero  level  surface  would  be  possible  using  equip¬ 
ment  of  this  nature  and  would  produce  results  similar  to 
what  has  been  obtained  but  in  a  more  graphic  and  versatile 


manner. 
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Even  this  relatively  simple  computer  assistance  in 
finding  spatial  solutions  to  the  given  existence  theorem 
demonstrates  the  possibility  of  using  computer  mapping 
programs  to  find  spatial  solutions  to  that  whole  class 
of  existence  theorems  having  explicit  or  implicit  spatial 
connotations.  After  initial  problems  related  to  adopting 
the  particular  existence  theorem  to  a  computer  solution 
are  resolved,  computer  mapping  techniques  are  capable  of 
quick  and  accurate  spatial  solutions  with  a  minimum  of 
human  manipulation.  Thi?  is  another  indication  that  the 
problem  solving  capabilities  of  computer  mapping  are  at 
least  as  diverse  and  effective  as  are  its  proven  capacities 
graphically  to  represent  stores  of  spatially  ordered 
information.  These  capabilities  l;c  largely  untapped, 
waiting  to  be  exploited. 
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